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■5W  quasisteady  model^dgrived  by  the  author  in  an  earlier 
paperdEJ^is  extended  to  fhclude  lateral  boundary  conditions. 

The  approach  is  to  first  assume  a  well-posed  complete  time- 
dependent  problem,  including  boundary  conditions;  it  is  stressed 
that  to  each  boundary  condition  there  should  be  associated  a 
precise  physical  assumption.  It  is  then  shown  that  the  quasi¬ 
steady  assumption  can  be  applied  consistently  to  both  the  partial 
differential  equations  and  the  boundary  conditions,  thereby  — .p  /via* 


FORM 
t  JAN  71 


EDITION  OF  I  NOV  «S  IS  OBSOLETE 


3ft, 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (*hen  Data 


VO  /  2  r  e  o 


I 


obtaining  a  wei-1  -osed  mathematical  model  with  time  scales  suitable 
for  large  scale  atmospheric  flows.  Three  types  of  conditons  are 
considered  at  the  lateral  boundaries:  1)  outflow,  2)  inflow  (driven 
velocity  is  specified  and  is  essentially  independent  of  the  internal 
flow,  3)  inflow  (passive)  -  inflow  is  created  primarily  by  the  inter 
flow  configuration .  The  upper  boundary  conditions  include  -the* 
two  derived^n — t’i'f',  the  continuous  and  discontinuous  boundary 
conditions,  and  a  third  condition,  which  is  designed  to  allow  the 
flow  to  propagate  independently  of  the  height  of  the  region. 

Numerical  solutions  are  obtained  for  various  test  cases,  and 
convergence  of  the  calculations  is  demonstrated.  One  sees  that,  although 
the  calculations  are  reasonable,  both  from  a  mathematical  and 
physical  standpoint,  the  various  calculations  differ  from  each 
other  very  significantly,  both  qualitatively  and  quantitatively. 

It  is  suggested,  therefore,  that  the  process  of  specifying  and 
evaluating  boundary  conditions  will  proceed  more  efficiently 
if  more  physical  understanding  is  obtained  in  regard  to  relatively 
simple  flows,  such  as  a  wave  entering  a  stationarv  flow  r.t  a  lateral 


[11  P.  Gordon,  "Quasisteady  Primitive  Equations  with  Associated 

Upper  Boundary  Conditions",  J.  Math.  Physics,  V.20,  April,  1979. 
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Lateral  Boundary  Conditions  for  Quasisteariy 
Atmospheric  Flows 


Section  I:  Introduction 

In  ref.  [l]  quasisteady  equations  were  derived  on  the  basis  of  scale 
assumptions  approximately  satisfied  in  large-scale  atmospheric  flow.  Several 
upper  boundary  conditions  were  also  derived,  but  the  lateral  boundaries  were 
severely  simplified.  The  underlying  philosophy  used  in  the  above  work  was  as 
follows:  Assuming  that  the  complete  time-dependent  problem  is  well-posed 
(this  includes  specification  of  boundary  and  initial  conditions),  derive  internal 
equations  and  boundary  conditions  bv  employing  clearly  stated  scale  assumptions. 

In  the  present  paper,  an  attempt  is  made  to  extend  this  approach  to  include  lateral 
boundary  conditions. 

Difficulties  with  boundary  conditions,  particularly  at  free  boundaries, 
arise  in  almost  all  areas  of  fluid  dynamics.  In  the  case  of  supersonic  flow, 
the  problems  are  generally  simpler,  since  one  can  expect  perturbations  at  such 
boundaries  not  to  feed  back  upstream  (see,  for  example,  ref. [2]).  Mathematically 
the  hydrodynamic  equations,  representing  a  hyperbolic  system,  are  reasonably 
well  understood.  This  is  particularly  true  for  the  one-dimensional  case,  although 
in  two  or  more  3pace  dimensions,  the  existence  and  unioueness  theorems  are  not 
so  well  formulated  [3].  It  is  in  fact  physical  understanding  that  is  lacking: 

It  is  perhaps  impossible  to  specify  physically  "accurate"  conditions  at  these 
free  boundaries.  This  situation. has  so  deteriorated  that  the  scientist  engaged 
in  computational  studies  feels  free  in  many  instances  to  ignore  the  physics  and 
mathematics.  Chen  [1*]  discusses  "...computational  boundary  conditions,  in  contrast 


to  physical  boundary  conditions  which  arc  originally  required..."  Also, 

in  his  survey  article  (5,  p.19],  Belotserkovski i  can  write  as  follows: 

"...the  basic  principle  underlying  the  statement 
of  the  conditions  is  that  no  substantial  disturbances 
should  penetrate  through  the  "open"  boundaries  of  the 
region  into  the  computational  region." 

Various  authors  have  written  of  the  necessity  of  neither  over-specifying 
nor  under-specifying  the  number  of  boundary  conditions  [6,7,8].  In  large 
scale  atmospheric  flow  an  additional  stumbling  block  presents  itself:  Because 
of  "simplifying"  assumptions,  such  as  the  hydrostatic  assumption,  the  equations 
are  no  longer  hyperbolic.  The  question  then  arises  anew  ns  to  even  the  proper 
number  of  boundary  conditions.  The  paner  of  Oliger  ft  fundstrom  [9]  would  seem 
to  indicate  that  there  is  no  "proper"  number:  A  major  result  in  this  paper  is 
that  the  hydrostatic  model  is  not  well-nosed  mathematically. 

The  above  result  has  been  quoted  often  in  the  last  few  years.  Unfortunately, 
at  least  in  the  author's  view,  the  result  has  at  times  been  interpreted  to  mean 
that  the  researcher  (since  the  equations  are  not  well-nosed  anyway)  can  now 
experiment  freelv  in  terms  of  boundary  conditions.  In  fact,  Oliger  and 
Hundstrom's  work  does  produce  the  expected  result:  Having  made  quasisteady 
assumptions  (the  hydrostatic  eouation  represents  such  an  assumption),  it  should 
not  be  surprising  that  not  all  boundary  conditions  are  still  compatible. 

This  turn  away  from  the  physics  is  disconcerting.  Historically,  mathematics 
has  been  productive  to  the  applied  community  when  closely  tied  to  the  physical 
problem.  Hadamard  [10],  who  seems  to  have  first  introduced  the  concept  of  a 
well-posed  problem,  emphasized  the  importance  of  maintaining  contact  with  the 
physical  problem.  In  [10,  p.3?l  he  writes  an  follows: 
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"...  But  it  is  remarkable,  on  the  other  hand, 

that  a  sure  guide  is  found  in  physical  interpretation: 

An  analytical  problem  always  being  correctly  set,  in 
our  use  of  the  phrase,  when  it  is  the  translation  of 
some  mechanical  or  physical  question..." 

Hadamard's  specific  concern  was  different  than  ours  (he  was  objecting  to 
the  mathematical  assumption  of  analytic  data),  but  his  point  is  still  valid. 

Tt  appears  to  the  author  that  in  the  area  of  computation  the  atmospheric 
physicist  ha3  relinquished  much  of  his  responsibility  to  the  mathematician/ 
numerical  analyst.  This  is  not  intended  to  disparage  the  latter  (if  anything, 
the  author  would  bo  so  classified).  Rather, the  intent  is  to  emphasize  the 
importance  of  at  least  partially  Judging  the  results  of  numerical  computations 
on  the  basis  of  the  Physical  problem  that  one  postulates  is  being  solved. 

For  example,  "reflective"  behavior  at  a  boundary  may  in  some  sense  be 
undesirable,  and  its  existence  may  in  some  sense  be  significant  for  diagnostic 
purposes.  However,  such  a  phenomenon  remains  at  best  a  synptom.  Since  one 
does  not  define  a  physical  problem  by  specifying  "non-roflect ive" ,  neither 
should  one  impose  such  a  criterion  on  the  mathematical  formulation.  And,  needless 
to  say,  neither  should  it  be  imposed  on  the  numerics. 

One  possible  procedure  for  specifying  outflow  boundary  conditions  is  to 
first  assume  physically  a  relatively  simple  form  for  the  solution  outside  the 
region,  and  then  specify  boundary  conditions  associated  with  this  flow.  (one 
boundary  condition  proposed  in  this  paper  is  of  this  form).  In  the  atmospheric 
case  such  a  procedure  is  difficult  because  so  few  "simple"  solutions  are  known. 
This  is  particularly  due  to  the  fact  that  in  an  atmospheric  environment  even 
simple  flows,  such  as  a  wave  entering  from  a  lateral  boundary  into  a  region  of 
stationary  flow,  have  extremely  complicated  solutions.  Nevertheless,  this  lack 


of  knowledge  is  highly  detrimental  to  both  the  process  of  specifying 
boundary  conditions  and  to  the  evaluation  of  corresponding  calculations. 

The  point  to  be  emphasized  is  that  the  boundary  conditions,  however 
specified,  must  be  interpreted  physically  in  terms  of  what  is  inherently 
being  assumed  regarding  the  external  flow  configuration.  This  philosophy 
does  differ  greatly  from  that  employed  by  most  atmospheric  scientists  engaged 
in  numerical  calculations.  Many  of  the  techniques  presently  in  use  (see, 
for  example,  Orlanski  [ll],  Perkey  &  Kreitzberg  [12],  and  Kemp  and  Lilley 
[13])  share  the  common  characteristics  o**  not  referring  to  the  external  flow. 

Boundary  conditions  are  discussed  for  the  time-dependent  case  in  Section  II 
and  for  the  quasisteady  case  in  Section  III.  Calculations  using  these  conditions 
are  described  in  Section  IV,  V  and  VI:  Section  IV  considers  a  problem 
involving  internally  generated  flow.  Section  V  considers  a  bell-shaped  wave 
entering  at  a  lateral  boundary,  and  Section  VI  considers  a  "flat"  wave 
entering  at  a  lateral  boundary.  Section  VII  summarizes  possible  conclusions 
of  the  study. 

An  accuracy  study,  in  terns  of  a  decreased  mesh  size,  was  conducted  for 
all  calculations.  Results  of  this  study  are  shown  for  the  bell-wave  calculation 


in  Section  V. 
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The  wUMitJ  to  h#  u*c\l  are  similar  to  those  defined  in  ref.  I,  except 
that  the  velocity  component  *  arc  nov  norma l  i  t  ed . 
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The  various  parameters  defined  above  are  as  follows: 

a  *  Cx/c.  *  r  ■  c  ♦  wc  ♦  uc  ,  y  =  1  ♦  R/c  .  c2  =  YRT, 
*  t  z  *  Vs 


^  *  r/p0*  ^  =  T/Tq»  ^1  =  position  of  the  lower 


boundary , 


f?  *  position  of  upper  boundary,  p  =  density,  T  =  temperature, 
w  »  vertical  velocity,  u  *  horizontal  velocity,  P  =  vas  constant, 
cy  «  specific  heat,  x  =  horizontal  distance,  z  *  vertical  distance, 
t  *  time,  g  *  acceleration  due  to  gravity,  p  *  pressure  *  pRT,  p^  and 


Tq  are  reference  values. 

The  eauations  for  the  floating  top  and  for  the  continuous  and  discontinuous 
boundary  conditions  (l:  eqs.  10,  15,  21]  take  respectively  the  following  fora: 


(f*).  -  w/I+F" 


(y-1) n  ~ 

IT  **  -  -»■  r -  t.Y 

T  C 

s 


(4.1) 


(A. 2) 


(  Y  -  1 )  n  ~  ~ 

w  :  w  <  0 


ir_- 1 )  *  s  . 

cs 


(4.3) 


:f2)T  :  w  >  0 


Remark:  In  ref.  1  the  discontinuous  boundary  condition  was  written  incorrectly 

_L 

in  that  the  scale  factor  (1  ♦  a2)  was  omitted  from  the  first  term  on  the  right. 


A  solution  is  to  be  obtained  in  the  region  0  <  t  <  1,  0  <_  n  <  1, 

0  <  r  <  t  , .  Before  considering  the  boundary  conditions  for  eq.(3)  at 
—  —  end 

the  four  spatial  boundaries,  a  general  condition  will  be  discussed. 

A  simplified  one-dimensional  model  consists  of  the  following  equations: 


c* 


"t  *  -  (y-1)*ux 


(5.2) 


Assuming  a  steady  wave,  moving  with  velocity  c^,  one  "watches"  the  wnve 
front  approach  bv  the  equation. 


dx 
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(5.3) 


Substituting  eq.(5.3).  cos.  (S.l)  and  (5.2)  take  the  following  form: 


e- 
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(5.4) 

(5.5) 


In  the  case  that  c  =  c  (that  is,  a  sound  wave),  cos.  (5.-*)  and  (S.S)  both 

w  s 

become  the  following: 


u  =  - - '- —  w 

<  (y-1)*  t  ‘ 


(S.p) 


Eq.(5*6)  has  been  used  as  a  boundary  condition  in  earlier  work  (lb,15l,  where 


the  following  physical  interoretat ion  was  given: 

At  a  boundary  eo.(^.f)  models  a  sound  wave, 
moving  into  an  infinite  region  of  undisturbed 
flow. 


Fq.(5*6)  will  be  used  in  this  paper  an  a  boundary  conditon.  An  obvious 
objection  would  be  that  the  problem  of  most  interest  in  atmospheric  situations 
will  not  be  sound-wave  dominated.  It  is  tempting,  therefore,  to  use  c^  in  eq.(S.3) 
"fudge-factor"  by  which  one  could  obtain  various  forms  of  either  eq.(5*^) 


as  a 


or  eq.(5.5).  This  approach  has  been  avoided,  and  in  this  regard  the  following 
points  are  pertinent: 

i)  The  philosophy  of  this  naoer,  as  discussed  in  the 
introduction,  is  that  the  boundary  conditions  must 
make  an  assumption  in  terms  of  the  external  flow.  It 
is  incorrect,  therefore,  to  calculate  c^  'adaptively'1 
from  conditions  of  the  internal  flow. 

ii)  If  specific  two-dimensional  steady-state  solutions 
to  atmospheric  problems  were  known,  then  rerhans  these 
solutions  could  be  used  to  derive  more  appropriate  boundary 
conditions. 

Specific  boundary  conditions  at  the  four  spatial  boundaries  will  now  be  given. 

A)  t  =  0  is  assumed  to  be  a  solid  boundary.  An  appropriate  boundary 
condition  is  w  =  0.  One  can  then  obtain  0  and  u  from  eqs.(3.1)  and  (3.2), 
while  it  can  be  calculated  from  the  following  eauation: 

Cs  c 

TTTJV’t  ‘N  "(  ‘  V  (7-0 

Eq .  (7*1)  represents  the  linearized  characteristic  variable  which  propagates 
information  from  £  >  0  to  C  =  0. 

B)  Because  of  ea.  (b.l),  £  -  1  is  a  streamline  and  conseouently  one  and 
only  one  boundary  condition  is  reouired  here.  The  three  differential  equations 
to  be  used  are  eqs.(3.l)  and  (3.2),  and  the  following  characteristic  equation; 

c  c 

WS  *■<  T?TJ?02  *  °3  >  •  (7-2) 

Three  upper  boundary  conditions  will  be  considered. 

l)  Eq.(U.2):  This  equation  was  derived  mathematically  in  ref  [1]  as  the 
continuous  boundary  condition.  The  interpretation  given  by  (6)  is  not  entirely 
valid  here  because  the  derivative  in  eq.(U.2)  is  with  respect  to  r  and  not  t. 
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Thus,  with  this  boundary  condition,  the  floating  tor  pushes  a  line  of 
constant  pressure  (  namely  r,  -  1  )  into  the  "undisturbed"  region  and  then 
superimposes  the  change  dictated  by  eq.(U.P).  In  essence,  the  characteristic 
variable  of  eq.  (7.1)  remains  constant  as  it  propagates  from  C  >  1  to  (  =  1; 
this,  of  course,  ignores  hydrostatic  changes  far  above  £  =  1 , 

2)  Eq.(U.3):  This  enuation  was  derived  in  [l]  as  the  discontinuous 
boundary  condition.  The  physical  interpretation  is  that  the  flow  at  c  =  1 
behaves  as  a  weak  compression  wave  moving  into  a  region  of  relatively  undisturbed 
flow.  However,  the  assumption  now  is  that  the  wave  acts  at  the  moving  boundary 
with  respect  to  the  undisturbed  hydrostatic  pressure  at  that  position.  That  is, 
the  change  dictated  by  eq.(U.?)  is  imposed  on  the  undisturbed  region  above 

C  ■  1;  the  second  term  on  the  right  in  eq.(U.3),  for  w  >  0,  accounts  for  the 
hydrostatic  change  in  pressure.  Interpretation  (6)  would  seem,  therefore,  to 
be  valid  for  this  boundary  condition. 

3)  Both  eqs.(U,2)  and  (b.3)  assume  that  the  flow  above  c  s  1  is  relatively 
undisturbed.  There  are  situations  where  this  physical  assumption  is  clearly  not 
justified.  One  such  example  would  be  that  of  an  incoming  lateral  wave,  with 
the  assumption  that  the  wave  is  also  entering  above  the  region  of  computation. 

A  possible  physical  assumption  for  this  flow  would  be  as  follows:  The  flow 
above  c  °  1  propagates  laterally  in  a  one-dimensional  fashion,  without  being 
significantly  affected  by  the  flow  below  5*1.  This  is  modeled  mathematically 
by  eq.(5.6),  with  interpretation  (6).  Transforming  to  t  ,  one  obtains  the  following 


x  + 
T 


(8) 


C)  At  the  lateral  boundaries,  three  boundary  conditions  are  required  for 


inflow  and  one  otherwise  (assuming  subsonic  flow).  At  a  non-inflow  point 


eq.(3.l)  can  be  used  along  with  the  following  linearized  characteristic 
equations: 

auT  -  wt  =  -(aG^-G^)  (9) 

S  1  \  *  -<«i  5  »  <“> 

In  eq.(io),  the  plus  sign  is  used  at  n  =  1  and  the  negative  sign  at  n  =  0. 

As  at  the  other  boundaries,  the  additional  boundary  conditions  must  reflect  a 
physical  assumption  regarding  the  outside  flow.  At  the  non-inflow  points  it 
was  decided  to  use  either  eq.(8),  with  the  same  physical  interpretation,  or 
the  condition  u  =  0. 

At  inflow  points,  eo.(lO)  -  with  the  appropriate  sign  -  can  be  used,  but 
three  additional  boundary  conditions  need  to  be  specified,  '"his  is  particularly 
troublesome  for  0  because  of  the  fact  that  perturbations  in  0  propagate  so 
slowly,  namely  at  the  flow  velocity:  th^  computational  times  to  be  considered 
in  this  paper  are  far  less  than  the  approximately  ll:  hours  required  to  traverse 
500  km.  at  lOm/sec.  It  was  decided,  therefore,  to  assume  that  0  is  not  perturbed 
by  outside  flow.  Physically,  this  would  imply  that  0  is  constant  on  incoming 
streamlines,  but  it  is  not  so  clear  how  best  to  interpret  this  mathematically. 

The  following  equation  was  eventually  settled  upon: 

Op  +  r  0  =  0 .  (11) 

In  order  to  determine  the  remaining  conditions,  it  is  necessary  to  distinguish 
between  "forced"  inflow  and  "passive"  inflow.  Forced  inflow,  which  represents 
physically  an  external  flow  being  impressed  on  the  region,  can  be  modeled  simply: 
Kither  velocity  or  pressure  is  specified.  For  purposes  of  the  present  paper,  it 


was  decided  to  specify  velocity: 
u  =  h^t.z) 

w  =  h2(t,z) 


(12.1) 

(12.2) 
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(Kqs.  (10)  and  (11),  t  hr»n ,  essentially  determine  0  ami  * ) ,  However,  the 
nee  ond  case,  where  flow  in  in  effect  pull  ml  into  the  region  because  of 
internally  produce*!  gradients,  in  more  complicated.  0m’  losallile  physical 
assumption  would  he  that  the  wave  enters  as  a  sound  wave  normal  to  the  surface. 
Mathematically,  this  translates  to  eqs.(Jl)  and  w  *  0, 

Summarising,  eq.(3)  ai  interior  points  and  the  rol loving  equations  at 
boundary  points  are  assumed  to  define  a  well-nosed  mathematical  formulation 
of  a  well-defined  physical  problem: 


1 ) 

C  ’ 

*  0: 

w 

-  0 ,  (1 

.1  ), 

(■!.?),  (7.1). 

(n.i) 

il) 

C  5 

•  It 

(■ l 

.1).  (1 

.?), 

(7.?)  and  one 

Of  ((>., 

.').(>». o.um.  (11.2) 

Hi) 

n  1 

•  0 

or  1 

• 

a) 

non 

-inflow:  ( 

1.1) 

.  (°  ) . ( 10) .  u  = 

0 

(l  i.  o 

or , 

(1. 

l).(o). 

(10) 

.(«). 

(n.»0 

h) 

for 

ceil 

i  11  f  1  OW  : 

(10),(U),(1?.1) 

.(!?.?) 

( it.'>) 

c ) 

mss  i  ve 

1  nflow 

:  (10), (11 ).(«).  w 

«  o 

(13.6) 

The  above  formulation  does  not  consider  the  corner  points:  (n  "  0,1;C  “  0,1). 
Several  severe  problems,  which  have  been  consistently  ignored  in  1  lie  literature, 
arise  at  these  points: 

t)  There  are  two  sets  of  equations,  corresponding  to  t  lie  two 
boundaries,  that  can  be  vised.  Generally,  one  attempts  t  o  use 
those  relat  ionships  which  are  forcing  the  '"low  configuration. 

For  example,  an  algebraic  boundary  condition,  because  of  continuity 
requirements,  in  considered  forcing. 

li)  01  nee  the  linearised  charset eri st in  equations  are  In  general  not 

valid,  1 t.  may  become  necessary  to  use  a  partial  differential  equation 

directly.  For  example,  if  the  flov  is  being  forced  at  n  •  0,  t  •  0, 

then  it  is  reasonable  to  impose  eqs.(l?.l),  and  (10.?),  with  ho?0, 

and  to  impose  0,  say  with  rq.(ll),  but  one  expects  to  calculate  * 

from  the  Internal  flow.  However,  neither  eqs,(3.?)  and  (7.1),  from 

the  t  *  0  set,  nor  eqs.(lO)  and  (II )  from  the  n  ■  0  set,  can  be 

applied.  In  this  ease  one  might  be  forced  to  use  eq.(l.t),  *  " 

t  • 

directly. 
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Section  III:  Specification  of  Boundary  Conditions  for  the  Qvm  si  steady  Model 


The  quasisteady  equations  are  obtained  from  eq.(l)  by  assuming  that  * 
and  w  are  in  quasisteady  equilibrium  with  respect  to  B  and  uf 3 ) . 

The  following  equations  result: 


rO  -  r  6  . 

C  L  n  ’ 


unchanged 


(lU.l) 


(]**.?) 


°2  *  °’ 


(l**.l) 


0^  -  0. 


(U.  *0 


It  remains  now  to  study  the  effect  of  the  quasisteady  assumption  on  the  boundary 

conditions.  The  underlying  principle  is  as  follows: 

If  a  boundary  equation  involves  or  w^ .  and  if  the  (15) 
equation  Involves  flow  conditions  internal  to  the  region 
of  computation,  then  the  conation  will  be  put  in  quasisteady 
equ i 1 i br i um . 

V 

At  C  =  0  the  situation  is  relatively  simnle.  w  *  0  and  rqs.(l.l)  and  (1.2) 
remain  unaffected  by  the  quasisteady  assumption.  According  to  (15),  however, 
eq,(7.l)  becomes  a  quasisteady  equation.  This  means  that  the  ripiht  side  of 
eq.(T.l)  is  set  to  zero.  Because  of  eqs.(l*i.l)  and  (l*t.!<),  then,  eq.(T.l)  is 


deleted  from  the  system. 

At  C  ■  1  a  similar  analysis  bolds:  fiqs.(l.l)  and  (1.8)  are  still  valid,  while 
eq.(7.2)  is  deleted.  On  the  other  hand,  the  upper  boundary  conditions  are  unaffected 
by  the  quasisteady  assumption:  Although  involving  quasisteady  variables,  these 
equations  do  not  require  internal  conditions  of  the  flow;  in  fact,  the  "derivation" 
of  these  equations  related  to  external  flow  conditions. 


At  the  lateral  boundaries  a  more  complicated  situation  exists.  First  note 
that  eqs.(lU.3)  and  (lU.U)  were  derived  only  at  internal  points  of  the  flow  and 
need  not  in  principal  he  imposed  at  the  boundaries.  However,  if  these  equations 
are  not  used  at  the  lateral  boundaries,  the  resulting  solution  may  exhibit  3teep 
gradients.  Such  gradients  will  be  inconsistent  with  the  scale  assumptions  upon 
which  the  quasisteady  equations  are  based  [1,  assumption  2],  (This  perhaps  relates 
to  the  discussion  in  the  introduction  reRardinR  the  results  of  Oliver  and  Bundstrom). 
One  concludes  as  follows: 

The  basic  quasisteady  assumption  fl,  assumption  0]  requires 

that  eqs.(llt.l)  and  (llj.lj)  bo  applied  at  all  Points  of  the  (l6) 

lateral  boundaries  except  Possible  at  the  corner  points. 

As  noted  in  ref.  [l],  the  effect  of  the  quasisteady  equations,  eqs.(]U.3) 
and  (lU.U),  is  to  remove  internal  time  dependence  from  the  variables  it  and  w. 

Boundary  conditions  must  account  for  the  time-dependence:  w  =  0  at  C  =  0 
accomplishes  this  for  w  and  the  upper  boundary  condition  provides  the  time 
dependence  for  « .  It  is  important  to  note  that  eqs.(U.P)  and  ( U . 3 )  provide 
time  dependence  relative  to  the  upper  flow,  while  eo.(B)  also  relates  to  the 
external  lateral  flow;  therefore,  if  eq.(ft)  is  used  at  the  top  it  is  necessary 
to  consider  the  possibility  of  incurring  inconsistencies  at  the  lateral  boundaries. 

Consider  first  non-inflow  points  and  assume  the  upper  boundary  condition  is 
either  eq.(b.2)  or  ( U .  3 ) .  If  the  lateral  equation  set  is  Riven  by  (13-3),  one 
would  use  eqs.(3.l)  and  u  =  0  for  the  quasisteady  model,  and  if  the  lateral 
equation  3et  is  Riven  by  (13. U),  one  would  use  eqs.(3.l)  and  (ft).  If  the  upper 
boundary  condition  is  Riven  by  eq.(ft),  the  situation  remains  the  same  for  equation 
set  (13.3)*  but  an  inconsistency  c.an  arise  with  respect  to  equation  set  (13. 1<). 
clearly  near  c  *  1,  eq.(fl)  cannot  be  used  to  calculate  u  "laterally”  and  n 
"vertically".  In  this  case,  since  it  seems  reasonable  to  assume  that  u  can  be 
obtained  from  the  internal  flow,  eq.(3.2)  has  been  used  directly  (this  assumes 
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that  the  required  time-dependence  at  the  lateral  boundary  is  provided  by  the 
upper  boundary  condition  in  conjunction  with  the  quasisteady  equations). 

At  forced  inflow  the  situation  is  reasonably  clear.  Eq.(lO)  is  deleted 
from  equation  set  (13. 5)  because  it  involves  internal  time  derivatives  of 
quasisteady  variables,  and  the  algebraic  equation  (12.2)  is  deleted  because 
it  is  inconsistent  with  the  quasisteady  equation,  eq.(ll<.3).  The  remaining 
equations,  eqs.(ll)  and  (12.2),  are  used  with  the  quasisteady  model. 

At  passive  inflow  points  eq.(ll)  can  still  be  applied  in  the  quasisteady 
case.  As  above,  eqs.(lO)  and  w  =  0  are  deleted.  The  remaining  equation  of 
equation  set  (1^.6),  namely  eq.(8),  can  be  used  if  the  unoer  boundary  condition 
is  given  by  eq.(U.P)  or  (b.3).  If  the  unper  boundary  condition  is  given  by 
eq.(8),  then  the  same  inconsistency,  as  discussed  above,  can  arise  if  eq.(8)  is 
also  used  at  the  lateral  boundary.  This  situation  has  not  yet  been  resolved. 
One  possible  solution  might  be  to  again  use  ea.(3.2)  for  u  (perhaps  it  is  in 
some  sense  the  "residue"  of  the  quasisteady  limit  process  relative  to  eq.(lOj. 
Further  study  is  required  for  this  case. 

Summarizing,  the  quasisteady  model  consists  of  eqs.(lb.3)  and  (lb.U)  at 
all  points  and  eqs.  (lU.l)  and  (lU.2)  at  all  points  except  possibly  n  =  0  or  1. 
w  =  0  is  imposed  at  c  =  0.  The  additional  feasible  boundary  conditions  to  be 
used  in  conjunction  with  these  equations  are  summarized  in  the  following  table: 


Top  Boundary 
Condition 

non- inflow 

|  forced  inflow 

passive  inflow 

(U.2)  or  (U.3) 

(3 .1 )  ,u  =  0 

or 

(3.1) ,  (8) 

(11), (12.1) 

1 

(11), (8) 

1 

(8) 

(3.1)  ,  u  =  0 

or 

(3.1) , (3. 2) 

(11), (12.1) 

! 

(11),  ? 

TABLE  1 


Section  IV:  Internally  Generated  Outflow  and  Inflow 

The  basic  numerical  scheme  can  be  found  in  ref.[l]  and  will  not  be 

repeated  here.  All  boundary  conditions  are  solved  implicitly  (  that  is,  at 

the  forward  time  step  ).  For  example,  at  a  lateral  boundary  eq.  (3.2)  would 

be  differenced  as,  — j—  (u?+y  -  u?  , )  =  )?+l  where  one-sided  differences 

At  \  l  % j  i  *«'/  t  i  f  j  « 

would  be  used  as  required.  Since  eqs.(lb.3)  and  ( 1 U . U )  are  also  solved 
implicitly,  the  iteration  problem  becomes  complicated. 

The  initial  data  used  for  the  test  problems  in  ref.fl)  consisted  of  a 
stationary  flow  with  a  lateral  pressure  gradient  in  a  region  o'*  height  h  and 
length  L.  Because  the  lateral  boundaries  were  assumed  to  be  solid  walls  (u  =  0 
was  the  boundary  condition),  the  flow  thereby  venerated  proceeded  to  reflect  from 
the  boundaries.  In  the  present  study,  the  right  boundary  will  be  a  free  boundary. 
The  lateral  boundary  conditions  must  therefore  allow  the  flow  to  move  through 
and  then,  as  it  turns  out,  to  move  back  into  the  region.  These  results  will  be 
compared  with  the  following  test  case:  The  region  is  tripled  to  length  3L, 
initial  data  in  0  <  x  <_  L  is  as  above,  data  for  the  additional  region  is  stationary 
flow  with  no  pressure  gradient.  This  latter  nroblem  is  solved  with  the  quasisteady 
model  and  u  =  0  at  the  lateral  boundaries.  During  a  reasonably  long  time  period, 
during  which  the  results  at  x  *  L  are  independent  of  the  boundary  conditions 
imposed  at  x  *  3L,  the  results  obtained  with  the  large  region  can  be  compared  with 
the  free  boundary  case. 

The  initial  flow  field,  for  the  region  o  x  <  1j8o  km.  *  L,  is  as  shown  in 
fig.  1  (this  is  reproduced  from  ref.fl]).  w  and  u  are  zero  throughout  at  time 
zero.  Overpressure,  with  no  pressure  gradient  in  l  <  x  <  3L,  is  ahown  in  fig.  2. 
For  this  teat  case  the  continuous  upper  boundary  condition,  eq.(l».2)  was  used  at 
C  ■  1,  and  u  ■  0  was  imposed  at  n  =  0  and  1.  Figure  3  shows  pressure  at  i  ■  0 
for  various  times.  Initially  a  "secondary"  pressure  front  forms  near  x  ■  L  (this 
can  be  seen  at  t  ■  1600),  but  by  ?U00  sec.  this  front  is  moving  toward  the  right 
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boundary.  This  pressure  distribution  is  accompanied  by  a  reverse  flow;  this 
can  be  seen  from  figures  U  and  5,  which  display  the  u  distribution  »t  C  *  0 
and  c  a  1. 

Ideally,  one  would  now  like  to  specify  boundary  conditions  at  x  *  L  so  that 
the  solution  thereby  obtained  would  reproduce  the  above  flow  in  0  <  x  <  L.  It  seems 
fairly  clear  that  for  a  problem  of  this  complexity  such  a  goal  cannot  be 
attained.  Nevertheless,  because  the  flow  is  internally  generated,  one  can 
expect  to  obtain  the  following  with  proper  boundary  conditions: 

a)  The  flow  should  be  qualitatively  correct  (in  this 
case  outflow  and  reverse  flow  should  occur  in  a 
generally  appropriate  time  scale). 

b)  Quantitatively,  the  results  should  be  comparable. 

If  inflow  occurs,  it  will  be  Passive  inflow.  Thus,  the  following  boundary 
conditions  are  used  (see  table  l):  Atx=n=0,  u=0  and  eq.(l.l),  and  at 
x/L  =  n  *  1,  eqs.(8)  and  (ll). 

This  solution  was  run  until  time  1*000  sec.  In  ref.  [l]  it  was  shown  that 

the  oscillation  time  for  this  boundary  condition  was  approximately  1600  sec. 

This  implies  a  wave  velocity  of  L/1600  ~  300  -?1—  ,  or  about  6^00  sec.  would  be 

sec. 

required  for  the  wave  to  travel  from  I.  to  3T>  and  back.  It  was  felt,  therefore, 
that  for  thi3  time  period  the  results  should  he  independent  of  the  boundary 
condition  at  x  *  3L. 

boundary  condition  (8)  is  intended  to  represent  a  situation  in  which  the 
external  flow  is  undisturbed.  Clearly,  for  this  problem  this  assumption  is  not 
valid  near  x  ■  L.  However,  it  perhaps  is  satisfied  for  the  flow  sufficiently 
distant  from  x  ■  L. 

Figure  6  compares  the  pressure  distribution  at  two  times.  The  solid  dots  are 
the  results  using  eq.(8)  at  x  ■  L.  At  time  800  sec.  the  comparison  is  very  good, 
while  at  time  U000  sec.  some  discrepancy  can  be  seen.  Figure  7  shows  the  comparison 


in  u  at  the  top  surface.  The  effect  of  eq.(8)  seems  to  he  to  delay  somewhat 
the  development  of  reverse  flow.  The  author  ,1udges  this  comparison  to  he  good: 

In  addition  to  maintaining  the  time  scale  of  the  problem,  the  results  remain 
reasonable  quantitatively  also:  It  might  perhaps  also  he  noted  that  this  is 
considered  to  he  a  difficult  test  case,  particularly  because  of  the  combined 
inflow-outflow  configuration  at  x  ■  L. 

It  seems  reasonable  to  suppose  that  even  better  agreement  could  be  obtained 

by  appropriately  "tuning"  the  factor  c  in  eq.(8).  For  the  reasons  discussed 

s 

in  Section  II,  such  a  procedure  was  not  at  tempted . 

A  calculation  was  also  made  with  eq.(3.?)  replacing  the  boundary  condition 
eq .  (8 ) ;  that,  is,  u  is  calculated  at  the  boundary  x  *  L.  The  Physical  significance 
is  not  clear,  particularly  because  this  condition  seems  to  make  no  direct  assumption 
regarding  the  external  flow.  It  is  interesting  to  note  t herefore  that  the  results, 
not  shown,  are  sensible  numerically,  but  3how  little  agreement  with  the  results 
obtained  with  the  large  region.  For  example,  at  time  h OOt'  sec.,  the  pressure 
distribution  has  oscillated  back  essentially  to  the  distribution  existing  at  time 


800  sec . 


Section  V:  Bell-Wave  Input  at  Left  Boundar 


For  this  case  the  initial  flow  is  stationary  and  hydrostatic:  u;v?0,  » 

3i 

^  ,  6(0, x,z)  *  6(0, 0,s)  and  6(0, 0,0  is  as  in  the  test  case  of  Section  IV 

o 


(l,  eq.(22.8)].  At  x  *  0  u  is  specified  as  a  function  of  time: 


f16  -wktc  -  -?y 


u(t,0,c)  =a 


0  <  T  <  C, 


(17) 


i^O:  *  >  V 

u  thereby  achieves  a  maximum  of  u  at  t  =  c,/2. 

max  1 

The  first  calculation  uses  eq.(3)  at  the  upper  boundary.  Consequently,  from 
table  1,  eq.(ll)  and  (12.1)  were  imposed  at  p  =  0  and  eos.(3.l)  and  (3.?)  at 
n  *  1  (note  that  inflow  at  n  s  1  does  not  occur  for  this  case).  The  values  of 
u  and  c.,  were  taken  respectively  as  2.5  and  3200.  Figure  3  shows  the  velocity 
distribution  at  the  top  surface  for  various  times.  At  time  U30C  sec.  this  wave 

has  left  the  region,  and  the  flow,  except  for  a  small  perturbntion  of  approximately 
has  returned  to  a  stationary  flow. 

This  residual  flow,  which  apparently  proceeds  to  oscillate  back  through  the 
region  with  values  |u|  <  .1  and  |w|  <  .005,  does  not  seem  to  be  numerical  error. 
Rather,  it  appears  to  be  related  to  the  difficulty  noted  in  Section  II  of  achieving 
a  true  steady-state  solution.  This  difficulty  is  due  to  the  slow  time  scale  of  6. 
At  time  6000  sec.,  the  variation  of  6  from  its  values  at  time  zero  is  still  on  the 
order  of  10  ;  this  is  sufficient  to  maintain  pressure  gradients  to  account  for 
the  values  noted  above.  To  partially  check  this  statement,  the  calculation  was 

repeated  with  6  held  constant  throughout  the  region  (6  .  6  =1.05).  The  "residual" 

o 

values  of  velocity  in  this  case  were  |u|  <  .001  and  [w|  <  .0001. 

The  shape  of  the  wave  form,  as  a  function  of  time,  is  shown  in  Figure  9.  The 
curve  labeled  x  *  0  represents  eq.(lT)  for  the  conditions  given  above.  The  curve 
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labeled  x  *  L  represents  the  velocity  curve  as  it  exits  the  region  at  the 


upper  boundary.  Since  the  region  was  initially  in  equilibrium,  one  expects 
the  wave  to  move  through  relatively  unchanged.  The  curve  shown  for  eq.(8) 
indicates  that  this  is  the  case.  In  addition,  there  is  no  damping  of  the  wave. 

The  exiting  velocity  curve  at  i  *  0  is  similar  to  that  shown  in  Figure  0,  except 
that  the  maximum  value  is  2.288. 

This  problem  was  next  solved  with  the  other  two  upper  boundary  conditions. 

From  table  1,  eqs.(ll)  and  (12.1)  are  used  at  n  *  0  and  eos.(l.l)  and  (8)  are 

used  at  n  “  1.  Again,  the  shape  of  the  curve  is  maintained  well.  The  vertical 

distribution  is  somewhat  more  pronounced  for  this  case.  For  example,  the  maximum 
outflow  velocity  at  t  *  0  is  u  =  2.151.  while  the  maximum  value  of  u  at  ;  *  1  is 

2.505.  Figure  (10)  shows  the  horizontal  velocity  distribution  at  c  *  0  for 

various  times. 

For  the  discontinuous  boundary  condition,  eq.(l».l',  the  results  are  markedly 
different.  Since  a  small  pressure  gradient  is  maintained  at  the  top  surface, 

horizontal  velocity  falls  off  rapidly  at  the  top.  In  fact,  as  noted  in  Figure  (11), 

the  wave  does  not  reach  the  far  boundary,  but  instead  decays  rapidly.  This  would 
indicate  that  a  large  amount  of  energy  is  being  dissipated  to  the  outside  urper 
region.  Velocity  at  c  *  0  does  build  ur  more  completely:  The  maximum  value  of 

u  at  (t  3  0,  n  *  *»)  is  1.5  and  occurs  at  time  2800  sec.,  while  the  maximum  value 

at  the  right  boundary  (1;*  0,n*l)  is  .1  and  occurs  at  time  l800  sec.  Fy  time 
U800  sec., a  complete  reverse  flow  exists. 

It  was  also  thought  of  interert  to  combine  the  upper  boundary  conditions. 
Eqs.(8)  and  (U.3)  together  make  little  sense,  since  each  contains  a  hydrostatic 
pressure  term.  However,  eqs.(8)  and  (b.2)  produce  the  ''ollowing  equation: 

.  .  btili  v  ♦  bell*  (S 

T  CT  C  T  T  C  T  C  V  to  J 

S  S 

Although  the  precise  physical  interpretation  of  rc.(l8>  is  not  completely  clear. 
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the  Idea  Is  to  combine  both  a  vertical  effect,  namely  eq. (»».,'),  with  the  horitontal 
effect  given  by  eo.(8).  Figure  12  shows  results  of  the  calculation  using  eq.(l8) 
at  the  top  and  eqs.(S.l)  ami  (3.?)  at  the  right.  lateral  boundary.  It  is  seen  that 
1-he  wave  form  is  much  like  that  obtained  with  the  individual  ecuat ions  (compare 
with  Figures  8  and  a).  However,  there  is  an  important  iistinction  in  that  eq.(l8) 
produces  some  damping  of  the  wave  form:  The  maximum  value  or  u  at  n  ■  1  is  2.  in 
and  occurs  at  approximately  t  “  1100  sec.  This  daraninc  appears  to  be  totally  an 
upper  boundary  effect,  unrelated  to  the  lateral  boundaries:  \’ote  that  the  wave 
is  already  depressed  by  time  2U00  sec.  Also,  comparison  at  an  earlier  time,  before 
the  wave  has  propagated  into  the  region,  shows  no  appreciable  difference  between 
the  three  cases  eq.(8),  eq.(4.2),  and  eq.(l8). 

The  three  boundary  conditions,  eqs . ( U .2) , ( 4 . 1) ,  and  (?',  differ  in  another 
imrortant  characteristic,  namely  dependence  on  height.  Vo. (8)  was  ’’designed"  to 
minimite  this  dependence:  The  physical  assumption  that  the  flow  is  also  being 
imposed  above  t  •  1  indicates  that  the  vave  profile  at  :  ■  1  should  be  relatively 
unaffected  by  the  height  of  the  lower  region.  This  was  verified  by  a  calculation 
with  a  height  of  5  km.  instead  of  10  km.  The  other  two  boundary  condition,  however, 
assume  undisturbed  flow  above  r.  *  1  and  consequently  should  be  height  dependent. 

This  dependence  is  reflected  in  the  form  or  the  equations:  Pis.  (!*..')  and  (4.3) 
involve  vertical  velocity,  which  is  a  quantity  very  sensitive  to  height. 

Figure  13  compares  two  calculations,  with  heights  of  10  km.  and  5  km. 
respectively,  both  using  eo.(l».2)  at  the  tee  and  the  same  lateral  boundary 
conditions  as  discussed  earlier.  Two  features  of  the  ''low  are  particularly 
interesting: 

1}  The  smaller  region  propagates  the  wave  more  slowly. 

This  is  perhaps  as  expected,  since  vertical  velocity 
is  smaller  and  consequently  pressure  builds  up  more 
slowly. 


?)  There  i s  a  damping  effect  related  to  height.  The 
5  Km.  case  produced  a  maximum  outlet  velocity  of 
1.89  m/  sec.,  while  for  the  10  km.  case  the  maximum 
outlet  velocity  was  shown  to  be  2.505.  Perhaps  for 
this  combination  of  wave  input  and  boundary  condition 
the  height  of  10  km.  is  some  kind  of  "magic  number": 

At  this  height  the  wave  propagates  across  the  field 
in  a  relatively  unperturbed  manner. 

Clearly,  there  is  a  significant  mathematical  relationship  between  the  wave 
input,  the  height  of  the  region,  and  the  upper  boundary  condition.  At  the 
present  time  what  one  requires  here  is  a  physical  explanation  of  this  relationship. 

As  discussed  in  ref.[l],  the  floating  top  is  an  important  part  of  the  model. 
Figure  lU  displays  the  motion  of  the  top  for  the  case  with  eq.(8)  as  the  upper 
boundary  condition.  In  general  this  motion  parallels  the  velocity  profile,  as 
can  be  seen  by  comparing  with  Figure  8.  This  was  also  true  for  the  cases  of 
boundary  conditions  (U.2)  and  (18).  The  maximum  perturbation  for  each  of  the 
three  cases  (using  eq . (8) , ( U .?) ,  and  (18)  respectively)  was  87.6  m.,  81.75m. , 
and  87. l*m.  For  the  discontinuous  boundary  condition,  eq.(l*.3),  the  top  boundary 
rose  21*1.5  m. 

As  in  the  test  case  of  Section  IV,  the  lateral  boundary  conditions  were 
tested  by  increasing  the  length  of  the  region.  The  region  was  doubled  in  length, 
with  the  boundary  conditions  formerly  imposed  at  x  =  L  now  being  imposed  at 
x  =  2L.  Figures  15  and  16  show  typical  results  of  this  investigation.  This 
particular  data  is  from  the  case  with  upper  boundary  condition  (18)  and  right 
lateral  conditions  given  by  eqs.(3.l)  and  (3.2).  The  solid  lines  are  the  results 
obtained  with  the  doubled  length,  while  the  solid  dots  are  the  results  with  the 
original  length  of  l*80km.  The  agreement  at  c  =  1  is  very  good,  while  some 
discrepancy  builds  up  a  5  =  0.  The  author  concludes  that  these  calculations 
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give  strong  assurance  that  the  lateral  boundary  conditions,  in  conjunction  with 
the  Riven  upper  boundary  condition,  are  consistently  modeling  the  hypothesized 
physical  situation.  In  this  case  the  hypothesized  physical  situation  is  that 
of  an  external  flow  which  is  undisturbed  "sufficiently  far"  from  the  region  of 
interest. 

Accuracy  of  the  computations,  relative  to  mesh  size,  is  considered  to  be 
an  important  part  of  a  numerical  study.  Of  the  three  classes  of  problems 
discussed  in  this  paper,  the  bell-wave  problem  of  this  section  displayed  the 
steepest  gradients  and  consequently  also  displayed  the  greatest  difficulty  in 
establishing  convergence. 

As  discussed  in  [l],  one  should  at  least  be  able  to  demonstrate  that  the 
solutions  of  a  numerical  algorithm  behave  as  though  they  are  part  of  a  convergent 
seciuence  of  calculations.  For  a  first  order  method  the  following  criterion  was 
derived  in  [ 1  ] : 


A2 

A1 


< 


•  5, 


(19) 


where  A1  =  f*(r,n,C.h)  -  ,  A?  =  f*(T,n,C.|-)  -  f*(r,n,C,£),  and 

f*(T,n,C,h)  is  the  numerical  approximation,  obtained  with  step-size  h,  to  the 
variable  f(T,n,c). 

For  the  problem  using  eq.(8)  at  the  upper  boundary,  the  following  sequence 
of  calculations  was  made: 

run  1:  Ax  =  20  km.,  Az~'1.25  km..  At  =  Uo  sec., 

run  2:  Ax  =  10  km.,  AZ  —  .625  km..  At  ■  20  sec., 

run  3:  Ax  »  5  km.,  Az--.  312*5  km..  At  «  10  sec. 

At  time  1600  sec.,  the  following  tables  display  typical  results  for  various 
quantities  in  the  flow  field: 
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In  general,  criterion  (10)  is  satisfied  very  veil.  Af^  represents  the 
actual  motion  of  the  top  surface  and  0  represents  the  deviation  of  6  from  its 
initial  value:  In  both  these  cases  there  is  some  question  regarding  the 
significance  of  the  digits  shown  in  the  table.  Note  also  that  P(x,0,l)  changes 
little  from  its  initial  value  (namely  ?6o.039),  while  P(r,0,0)  changes  much 
more  significantly  (from  its  initial  value  of  10?°. 792). 


Section  VI:  Flat  Wave  Input  at  Left  Boundary 


The  problem  in  precisely  as  that  discussed  in  Section  V,  except  that 

eq.(l7)»  the  input  function  for  u,  in  replaced  by  the  following: 

? 

u  ( T  ,  0 ,  C )  3  u  (l-e‘ClT  )  .  (?0) 

max 

For  the  sample  calculation  u  was  ngain  ? and  c,  *  10-<\  Using  the 

max  1 

upper  boundary  condition  eq.(fl)  and  the  lateral  boundary  conditions  (at  n  3  1 ) 
eqs.(3.l)  and  (3.S),  one  obtains  the  solution  depicted  in  Figures  17  and  18. 

T)\e  solution  at  the  too  surface  in  as  expected:  the  velocity  builds  up, 
according  to  eq.(20),  and  propagates  across  the  field.  However,  the  solution 
at  c  =  0  was  not  expected,  in  that  the  maximum  velocity  at  n  =  0  did  not 
propagate  in  the  expected  time  scale.  (By  time  6000  sec.,  the  profile  changes 
little  from  that  at  time  it 000  sec.).  The "discrepancy"  is  approximately  the 
same  as  that  seen  earlier  in  the  bell-wave  calculation. 

Tt  is  clear,  from  simple  conservation  of  mans  considerations,  that  the 
solution  as  shown  cannot  be  a  steady-state  solution.  It  is  also  clear,  from 
calculations  with  reduced  mesh  size,  that  numerical  error  cannot  account  for  the 
unexpected  form  of  the  solution.  One  next  needs  to  ask  whether  the  calculation 
satisfies  the  original  time-dependent  problem.  Thin  is  difficult  to  .judge,  hut 
if  the  answer  is  no,  then  presumably  the  quasisteady  equations  are  not  valid. 

As  noted  in  ref ( 1  1 ,  t.be  quasisteady  equations  do  introduee  errors  into 
the  solution.  Tn  the  present  problem  this  error  can  be  seen  clearly  at  n  3  0: 

The  velocity  input  at  n  3  0  would,  with  eos.(l),  produce 
a  pressure  increase  with  a  significant  non-hydrostat ic 
vertical  variation  (at  least  in  t.he  short  term).  One 
expects,  over  the  long  term,- that  the  vertically  propagating 
"fast"  characteristics  will  dissipate  the  "non-hydrostatic" 
variation.  The  quasi  steady  equations,  on  the  other  hand,  ignore 
the  short-term  time-development,  and  thereby  inherently  assume 


that  any  steep  gradients  at  n  ■  0  are  not  important 
in  the  long  term. 

It  is  not  thought  that.  the  above  error  would  account  for  the  result* 
shown  in  Figures  17  end  18.  One  expects  the  error  introduced  by  the  quasi  steady 
equations  to  be  mani Tented,  not  in  the  final  flow  nrofilea,  but  in  the  short¬ 
term  time  scale.  Helative  to  the  time-acale  of  the  problem,  this  error  should 
be  smal l . 

Referring  to  eq.(l.r),  one  sees  that  the  x  term  Is  the  most  important 

n 

factor  in  the  propagation  of  u.  The  coefficient  of  this  term  involves  only  0. 

A  calculation  with  0  constant  should,  then,  remove  much  of  the  vertical  variation, 
hut  would  not.  affect  the  auaslateady  assumptions .  Puch  a  calculation  was  made, 
with  0  5  l.Oh.  In  thin  case  the  variation  In  sound  arced  wan  even  more  pronounced 


(c  ) 

s  bottom 


7 — v -  -  l..'?.  The  velocity  profiles  ore  much  like  those  of  Figure  If, 

^  s  top 

except,  slightly  retarded.  However,  there  is  no  significant  variation  between 


the  top  and  bottom  horizontal  velocity.  Furthermore,  this  calculation  achieves 
a  steady-state:  by  time  UU00  sec.,  the  variation  in  (he  position  of  the  top 
surface  is  1.1m  (from  n  x  0  to  r\  ■  l) ,  and  the  maximum  value  of  |w|  ts  ,0001. 

A  more  careful  analysis  of  the  original  calculation  (with  variable  <0  shows 
that  the  solution,  although  not  in  steady-state,  is  very  slowly  varying.  For 
example,  at  time  6000  r.ec .  the  variation  in  the  position  of  the  top  surface  is 
5.5?m.,  and  the  maximum  value  of  |v|  is  ,00b?.  In  light  of  the  constant  0 
calculation,  it  would  appear  that  the  variables  u,v,  and  s  are  In  qussi-ateady 
equilibrium  with  respect  to  the  0  field.  (This  problem  of  the  slowly  varying  R 
wan  discussed  in  Pec t ion  IT). 

If  one  then  assumes  that  the  quasi  steady  model  is  producing  a  valid  solution 
to  the  time-dependent  problem,  one  then  needs  to  ask  whether  or  not  the  solution 


is  physically  meaningful .  As  dismissed  earlier,  this  question  has  not  yet  been 


answered.  If  the  answer  is  yen,  one  could  develop  reasonable  confidence  In 
the  boundary  conditions.  If  the  answer  is  no,  then  one  can  reexamine  various 
features  of  the  model,  such  an  the  velocity  input  function,  quasisteady  forms 
of  the  characteri Stic  equations,  and  boundary  conditions.  In  any  case,  physical 
significance  of  the  solution  needs  to  be  considered. 

This  problem  was  also  solved  with  boundary  conditions  ( U .  1' )  and  ( U . 3 )  at. 
the  top;  eq8.(3.l)  and  (fl)  were  used  at  the  right  boundary.  With  eq.(U,2)  the 
wave  propagates  across  the  region  in  a  reasonable  fashion,  but  a  significant 
overshoot  in  velocity  (approximately  10*)  occurs  at  the  right  corner  point.  It 
is  not  yet  clear  whether  the  problem  is  due  to  the  boundary  conditions  or  to  the 
corner  point  equations. 

With  eq.(h.3)  the  wave  propagates  across  the  field  at  a  much  slower  velocity. 
This  can  be  seen  from  Figure  19,  which  displays  the  velocity  profiles  at  c  *  0 
for  the  two  cases  at  a  fixed  time.  The  solution  with  eq.(U.l)  displays  another 
interesting  characteristic:  Rather  than  achieving  a  steady-state,  the  flow 
continues  to  push  upward  into  the  undisturbed  flow.  Figure  90  shows  the  horizontal 
velocity  at  C  3  1.  Near  n  =  0  the  flow  appears  to  achieve  a  linear  profile 
rather  than  a  constant  profile.  This,  then.  Produces  a  non-zero  vertical  velocity 
near  n  *  0,  and  this  in  turn  allows  the  flow  to  push  upward  into  the  region. 

This  upward  motion,  in  terms  of  the  position  of  the  tor  surface,  is  shown  in 
Figure  21.  It  is  conceivable  that  an  appropriate  incoming  lateral  flow  could 
move  vertically  into  the  less  dense  fluid,  rather  than  propagate  horizontally. 
Iloundary  condition  (I*. 3)  could  be  suitable  for  this  typo  of  flow. 


Section  VII:  Cone lun Ions 

For  the  most  part  this  study  is  based  on  the  following  two  assumptions: 

l)  A  boundary  condition,  by  its  very  nature,  mvist  reflect 
a  physical  assumption  in  regard  to  the  external  flow 
configuration.  Consequently ,  a  proposed  physical  interpretation 
should  be  associated  with  any  proposed  boundary  condition. 

P)  Fq.(3)  i 3  to  be  solved  in  a  region  of  height  h  and  length  L,  with 
h<<L.  It  is  assumed  that  boundary  conditions  and  initial 
conditions  to  be  imposed  on  Fq.(t)  are  such  t hat  the  flow 
variables  will  experience  significant  variations  only  over  time 
scales  which  are  large  compared  to  h/L. 

Remark:  The  second  of  the  above  is  the  basic  scale  assumption,  and  vas  formulated 
in  ref.ll]  as  assumption  (P). 

A  primary  purpose  of  this  paper  was  to  investigate  the  validity  of  the 
following  hypotheses: 

l)  There  exists  a  significant  lack  of  Physical  understanding  in 
regard  to  several  relatively  simple  flows  in  an  atmospheric 
environment.  This  situation  seriously  hampers  the  process  of 
specifying  and  evaluating  boundary  conditions  to  be  vised  in 
mathematical  models  for  large-scale  atmospheric  flow. 

ii)  Assume  that  Fq.(d),  together  with  a  given  set  of  boundary  and 
initial  conditions,  defines  a  well-posed  mathematical  problem. 

It  is  then  possible,  by  consistently  applying  the  above  assumption 
(P)  to  both  the  partial  differential  eouat. ions  and  the  boundary 
conditions,  to  obtain  well-posed  mathematical  models  whose  time- 
dependent  scale  is  suitable  for  large-scale  atmospheric  flow. 

In  support  of  the  first  hypothesis,  the  following  specific  points  are  noted : 

a)  because  of  the  slow  reaction  time  for  0,  it  appears  that 
actual  steady-state  solutions  cannot,  be  reached  in  the  time 
scale  of  these  problems.  This  behavior  of  0  needs  to  be 
factored  into  the  physical  interpretation  of  the  problem. 


b)  In  the  problem  of  a  wave  enteriag  at  a  lateral  boundary, 
several  factors  may  significantly  affect  the  profile  of 
the  propagating  waves.  These  factors  include  the  height 

of  the  region  of  entry  and  the  choice  o*  the  unner  boundary 
condition.  Mathematically,  all  solutions  seemed  reasonable. 
Consequently,  further  evaluation  requires  more  physical  tnpiit 
regarding  definition  of  the  problem. 

c)  In  Section  IV  an  example  was  described  in  which  the  solution, 
although  reasonable  mathematically,  did  not  agree  with  the 
solution  for  the  extended  region.  Again,  physical  input 
(which  was  not  postulated  for  this  particular  set  of  boundary 
conditions)  is  required. 

support  of  the  second  hypothesis,  the  following  nnecific  points  are 


a)  The  quasisteady  model  operates  on  a  time-scale  suitable  for 
large-scale  atmospheric  flow. 

b)  All  solutions  were  shown  to  be  stable  and  continuously 
dependent  on  the  data.  (Tn  ref.[l]  it  was  indicated  that 
convergence  in  terms  of  a  decreasing  mesh  size  is  a  test 
of  staoility). 

c)  All  solutions  shown  appear  to  have  a  sensible  physical 
interpretation.  However,  in  one  problem  (see  Cection  VI) 

an  overshoot  occurred  at  a  corner  point;  this  might  indicate 
that,  for  this  choice  of  boundary  conditions,  the  quasist.eady 
limit  was  not  taken  correctly. 

d)  The  proposed  lateral  boundary  conditions  behaved  well  and  it 
was  shown  that  they  could  be  interpreted  physically.  In 
particular,  the  comparison  with  solutions  for  an  extended 
region  was  very  satisfactory. 
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FIGURE  2:  OVERPRESSURE  AT  INITIAL  Tlltmr  -  p(o,x,t)  .  P(0,t,,) 
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FIGURE  5:  Vat  top  surface 
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FIGURE  17:  HORIZONTAL  VELOCITY  AT  TOP  SURFACE 
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